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The properties of a dilute electron gas, coupled to the lattice degrees of freedom, are studied and compared 
with the properties of an electron gas at half-filling, where spinless fermions with two orbitals per lattice 
site are considered. The simplest model which includes both the local electron-lattice interaction of the 
Jahn-Teller type and the electronic correlations is the E ® /3-Jahn-Teller-Hubbard model. We analyze 
the formation and stability of Jahn-Teller polarons and bipolarons, respectively. Our approach is based on a 
hopping expansion in the strong-coupling regime. The results are compared with recently published findings 
for the Hubbard-Holstein model 1 1 2 1. The special case of the Jahn-Teller-Hubbard model at half-filling is 
mapped on a spin- 1/2 Heisenberg model with phonon-dependent coupling constants. This has been derived 
within a projection formalism that provides a continued-fraction representation of the Green's function. We 
study the exact solution for two and three particles and compare it with the effective theory on the infinite 
lattice with one particle per site. 



1 The Model The Hamiltonian for fermions with spin a and pseudospin j — 9,e, coupled to phonons, 
is given hy H = Ht + Hq, where Ht is a nearest neighbor tunneling term for the fermions 

Ht^-t J2 J2 4i^^yy- + ^-c- 

<jJ'>'T=Ta7=e,c 

and Hq is local, containing a (Hubbard) interaction and a term for dispersionless phonons of energy luq 
with Hq = and 

^^oj = wo&jfcj + gib] + bj)'Yiniea ± n^ea) + XI ^"J'rT'^jTi + C^o(«jeT + f^jei)(?^jeT + ^j^i)' 

(7 7 

where the plus (minus) sign refers to Holstein (E ^ (3 Jahn-Teller) electron-phonon coupling 0. Hoj is 
diagonaUzed by a Lang-Firsov transformation |3 and has energies 

Eoi = ujouj - Epl^^in^g^ ± nj„) + ^ C/nj^T'^jTi + Uo{nie^ + ?^jei)('^j£T + ^iei) (1) 

(7 7 

if there are nj(> 0) phonons and njj„[= 0, 1) electrons with (pseudo)spin a (7) at site j. Each electron 
has an energy gain Ep = The regime U < 2Ep has an attractive interaction, leading to an on-site 

bipolaron. 
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2 The Bipolaron Problem To study the binding properties of polarons we consider two electrons on a 
lattice with two orbitals {6 and e) per lattice site. In a strong coupling approach we study the similarities 
and differences of the non-degenerate Holstein-Hubbard (n^ or ng — 0) (HH) and the E ® /3-Jahn-Teller- 
Hubbard model (JTH). We restrict ourselves to the case of repulsive on-site interaction (U > 2Ep). There- 
fore the ground state of Hq has single occupied sites. If both electrons are in the same orbital, the situation 
is the same as in the HH due to the orbital conserving hopping. We show in the following, however, 
that the situation for electrons in the 9 and the e orbitals is different. This is a consequence of the fact that 
we have a different distortion configuration associated with different orbitals. 

For a hypercubic lattice a hopping expansion leads to an exponentially reduced hopping rate te~^p/^° 
in first order. The results in second order for each possible hopping process can be written in general as 



The first series is the result of processes with empty and singly-occupied sites as excited states, while the 
second series takes doubly-occupied states into account. 7(0, b) is related to the incomplete Gamma func- 
tion fSl. For electrons in the same orbital a is equal to [U — 2£'p)/ajo and for different orbital occupation 
a is given as {Uq + 2i?p)/(jJo- Values of b are ±Ep/u!o and ±.2Ep/uJo. The positive sign corresponds 
only to processes that describe hopping in different orbitals. One process takes place in the 9, the other 
in the e orbital. Hence they only appear in the JTH but not in the HH. The terms related to the processes 
explained above are related to processes that change the orbital occupancy at least at one site. This means, 
for instance, that the initial state has an electron in orbital 9 and the final state an electron in orbital e at 
one site. Only the electrons at nearest neighbor sites are available for these processes. If both hopping 
processes take place in the same orbital the minus sign appears. These are the only contributions in the 
HH. They also exist in the JTH if either both electrons are in the same orbital or if only one of the two 
electrons moves. The factor 2 arises if only two sites are involved in the process. While the series with 
6 > reduces the exponential prefactor even further, the series with b < cancel it in part or, in the case 
of —2Ep/ojo, even completely. This means that hopping processes that do not change the lattice configura- 
tion, i. e. the orbital occupancy, are not suppressed exponentially in the strong coupling regime. Processes 
with exchange in the same orbital are included in this case but not those in different orbitals. The latter 
change the orbital occupancy and therefore the lattice configuration. In a strong coupling approach we 
neglect the exponentially reduced terms. For nearest-neighbor electrons in the same orbital we have the 
possibility of an exchange, which eventually leads to a nearest-neighbor singlet and triplet configuration as 
in the HH. The singlet, a nearest-neighbor bipolaron, is the ground state if J7 < AEp. The corresponding 
binding energy does not dependent on the coordination number of the lattice. This is an artefact of our 
approximation, and it is A = t'^/Ep — At^ /U . On the other hand, the exchange of electrons in different 
orbitals vanishes in the strong coupling regime. The lowest energy state of electrons in different orbitals 
is that of unbound polarons, and it is equal to that of unbound polarons in the same orbital. The nearest 
neighbor states are higher in energy with a finite energy shift with respect to the unbound states. 

3 Jahn-Teller Effect for Systems with One Particle per Site For simplicity, the spin of the electron 
is ignored now. Formally, the orbital has a similar meaning as the spin in the sense that it appears as an 
additional degree of freedom. This gives us orbital-depending effects in tunneling, very similar to spin- 
dependent effects 1 6 1 . In this case the ground state for sufficiently small hopping rate t compared with Ep, 
ujQ and Uo is a singly occupied lattice without phonons. 

In the following we shall consider the Green's function G{z) = (z — H)^^ and its projection to the 
subspace with one fermion per lattice site and no phonon Pq. The Green's function satisfies the operator 
identity |7| 



Po{z-H)-^Po = iPo{z-H)Po-PoHPiiz-H)^^PiHPQ)^^ = {z- PoHqPo- H^f f)o\ (3) 




and 




(2) 
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where Pi projects to the complement: Pi = 1 — Pq. {■■■)J^ is the inverse on the projected Hilbert 
space. Since the Hamiltonian H preserves the total number of fermions in the system, the matrix of 
H in fermion-number representation is given by diagonal blocks. Consequently, the projected matrix 
Hef f = PoHPi {z - H)];^PiHPo acts only on the Hilbert space with a total fermion number equal to the 
number of lattice sites. In other words, all virtual processes in H^f / are creation and annihilation processes 
of pairs of doubly occupied and empty sites. Equation Q is a recursion relation on the projected Hilbert 
spaces and can be iterated. It provides a continued-fraction representation Q. A truncation after the first 



iteration by replacing H 
model 



Hq in (z — H) ^ leads to the approximation of Hef f by the XXZ Heisenberg 



^e//«^ E [aniS!S^,-i) + an{Sfsp + sysy) 



<jj'> 

with coupling constants for a lattice with sites 



(4) 



a\i{z) = -4- 



-2EplujQ 



-1 



LUq 



-7 



Uo 



-{N 



UJQ 



and 



2^Ep 2Ep 



UJQ 



For weak and strong electron-phonon coupling we obtain the same behavior for a^^: 

4 

^ z+{N- 2)Ep - Uo 
but a different behavior for a^i, where 

an ^ '^ll (ff/^o < 1), an ~ (s/t^o > !)• 



(5) 



(6) 



(7) 



Thus the weak-coupling limit of Hef / gives an isotropic Heisenberg model, whereas the strong-coupling 
Umit leads to an Ising model. 



3.1 Three-atomic Molecules As a simple example we study the effect of the geometry on properties 
in the case of a molecule with three atoms (A^ = 3), either in a stretched (SM) or triangular (TM) configu- 
ration. The Hilbert space of the full Hamiltonian has infinite dimensions, whereas the Hilbert space of the 
Hamiltonian of the projected Green's function with one electron per site and no phononhas only dimension 
d = 8. The problem of diagonalizing the effective Hamiltonian is further reduced by the global spin-flip 
symmetry to d = 4. However, the eigenvalues of Hef / are complicated functions of z due to the virtual 
hopping processes which include the creation of an arbitrary number of phonons. The projected Green's 
function of Eq. (|3} can be expressed in terms of eigenstates and eigenvalues of the effective Hamiltonian 
Hef f as {Ej I (z — H)~^\Ej) = l/{z — Ej (z)). The energy functions Ej (z) for the stretched and for the 
triangular molecule are listed in Tabled The geometric degeneracy of the TM gives a doubly degenerate 
ground state for Eq. Therefore, the energy E2 is absent in this case. 



3.2 Discussion of the pole structure The energy levels of the molecules are poles of the projected 
Green's functions l/(z — Ej), i.e. they are solutions of the equation z = Ej{z). The parameter z appears 
in a|| and a^i (and therefore in Ej) only in the combination z' — z + Ep. In general, we expect a complex 
solution z' — X + iy for z — Ej. However, it turns out from the properties of the incomplete Gamma 
function that the imaginary part is always y = and that the real part satisfies an equation of the form 

f c 

x^—S^ — (cm > 0, am^i^Qm-2Ep + Uo)- (8) 

2 ^ X - am 

m 

Consequently, there are solutions Xm with < Xm < am + wq. 
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Table 1 Energies for the TM and the SM with Table 2 Values of the pa- 

Ao = -fait - iy^^Ts^^andAa = -fait + ^^J^~hS^^. rameter A in Eq. |9}. 



3.2.1 Pole structure in the asymptotics In weak coupling (WC) as well as strong coupling (SC), the 
Green's functions have only two poles for each Ej. These poles are solutions of a quadratic equation. To 
study all elements of the Green's functions in one case, the parameter A is introduced such that the Green's 
function reads 

{E,\{z - ny^lE,) = ^ , {z' = z + 3Ep) (9) 

^ z'-2Ep-Uo 

with the values of the parameter A given in Table|3 The poles in Eq. are 

z'1/2 = Ep + Uo/2 ± ^{Ep + Uol2f + 2f^\. (10) 

The ground state is related to m.in{z'^, zi^}, i.e., it is z!^. Moreover, A must be maximal. Thus Eq is the 
ground state, except for TM/SC, where the ground state has an additional degeneracy due to Eq — Ei. 

4 Conclusions It is argued in 1 2 1 that in the strong coupling regime the main source of bipolaron for- 
mation is the non-exponential off-diagonal matrix element in second order related to the exchange of 
neighboring electrons. In the E ® (3 case we found exponential decay for this exchange in the situation of 
different orbital occupancy. Therefore the unbound state is preferred. For one particle per site we showed 
that H^f f of the projected Green's function yields an XXZ Heisenberg model. Furthermore it should be 
noted that the corresponding Holstein-Hubbard model leads to an isotropic Heisenberg term. This was also 
discussed for the bipolaron problem 1 1 2 1. As an application we have discussed the spectral properties for 
stretched and triangular molecular configurations of three atoms. 
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